Algebraic Classification of Numerical Spacetimes and Black-Hole-Binary Remnants 



Manuela Campanelli, Carlos O. Lousto, and Yosef Zlochower 



o 
o 

>>■ 



1^ 



> 
O 

o 

m 



00 

o 



X 



Center for Computational Relativity and Cravitation, and 
School of Mathematical Sciences, Rochester Institute of Technology, 
78 Lomb Memorial Drive, Rochester, New York 14623 
(Dated: May 7, 2009) 

In this paper we develop a technique for determining the algebraic classification of a numerically 
generated spacetime, possibly resulting from a generic black-hole-binary merger, using the Newman- 
Penrose Weyl scalars. We demonstrate these techniques for a test case involving a close binary with 
arbitrarily oriented spins and unequal masses. We find that, post merger, the spacetime quickly 
approaches Petrov type II, and only approaches type D on much longer timescales. These techniques, 
in combination with techniques for evaluating acceleration and NUT parameters, allow us to begin 
to explore the validity of the "no-hair theorem" for generic merging-black-hole spacetimes. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70. Bw 



I. INTRODUCTION 

The recent breakthroughs in numerical relativity [11, [3, 
0] that allowed for stable evolutions of black-hole-binary 
spacetimes led to many advancements in our understand- 
ing of black- hole physics, and it is now possible to accu- 
rately simulate the merger process and examine its ef- 
fects in this hi ghly non-linear regime [1, H, |^, 0, H, H, 
[13, [m, [H [il IM d, M, M, m- Black-hole binaries 
radiate between 2% and 8% of their total mass and up 
to 40% of their angular momenta in the last few orbits, 
depending on the magnitude and direction of the spin 
components, during the merger [1,[1,[^ (ultra-relativistic 
head-on black-hole mergers can radiate up to ~ 14% of 
their mass [13 )• In addition, the radiation of net linear 
momentum by a black-hole binary leads to the recoil of 
the final remnant hole ^,'2T, 22, 23, 2^j2pj2ij23j2i, 

H, H [M 113, [H H [sa, 21, M, M, Efl ErM Ef , 

which can have astrophysi cally observable irnpqrtant ef- 
fects [H m, m, H SI, Si, S3 Si, Si, Hi, HH, M Hi and 

represents a possible strong-field test of General Relativ- 
ity (GR). 

In addition to important astrophysical applications, 
the two body problem in GR is intrinsically interest- 
ing because it provides the framework for analyzing the 
behavior of the theory in the highly-nonlinear, highly- 
dynamical, non-symmetrical regime. For example, the 
cosmic censorship hypothesis, that states that singular- 
ities in the universe should be cloaked by a horizon is 
under active investigation [1, [i, [i, [s^, [s^. In this pa- 
per we are interested in verifying the "no hair theorem" , 
which states that all black holes eventually relax into 
a state that can be described by three parameters, the 
mass, spin, and charge. Hence, the final merger rem- 
nants from multi-black-hole mergers [Hi, [s^l should be 
Kerr black hole [s^. 

The problem of determining the geometry of the final 
stage of a black-hole binary merger arises as a practi- 
cal question even in perturbative techniques, such as the 
Lazarus approach [53.[60j. which used a combined numer- 



ical and perturbative approach to simulate the waveforms 
from a binary merger. In the context of the Lazarus 
approach, it is crucial to determine when the transi- 
tion from numerical to perturbative evolutions is possi- 
ble, i.e. when the full numerical simulation could be ap- 
proximated by (relatively small) perturbations of a Kerr- 
rotating black hole, and a diagnostic, the S-invariant [6l| 



S = 27 J^/I^ 



(1) 



that is identically 1 for a Kerr spacetime, was developed 
to measure the closeness of the spacetime to an alge- 
braically special type II. However, the S-invariant by it- 
self is not sufficient to demonstrate that the spacetime is 
near Kerr because it does not distinguish between type 
II and type D spacetimes, nor does it imply that the 
acceleration and NUT parameters vanish. 

More recently, with the availability of new long term 
evolutions, one of the consistency tests performed is the 
agreement of the total angular momentum of the system 
when computed in three different ways: by measuring 
the angular momentum (and mass) of the remnant black 
hole [3 M, HI] using the isolated horizon formulae [6^ . 
by measuring the total energy and angular momentum 
radiated (63. , ,64] and subtracting it from the total initial 
values, and by looking at the quasi-normal frequencies of 
the late-time waveforms and associate them with those 
of a rotating Kerr hole with mass M and angular mo- 
mentum per mass a [li] . The rough agreement of those 
values represents indirect evidence that the final black 
hole is of the Kerr type. Furthermore, in Ref. [gi], where 
the authors of that paper presented very-high-accuracy 
waveforms from the merger of an equal-mass black-hole 
binary, it was shown that the minimum and maximum 
values of the scalar curvature on the remnant horizon 
agreed with the Kerr values. 

No hair theorems assume a stationary Killing vec- 
tor 58] as characterizations of the Kerr geometry [66l.f67t. 
While one can classify spacetimes based on their symme- 
try properties, here we will use a classification method 
based on the algebraic properties of generic spacetimes 
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without a-priori assumptions about symmetries. 

Demonstrating that the remnant of a black-hole 
merger approaches Kerr asymptotically (in time) would 
also help answer open questions about the stability of 
Kerr under arbitrary perturbations. The stability of the 
Kerr spacetime under linear perturbations has only been 
proven mode-by-mode [6^], and the interior of the hole 
may even be unstable [69i]. Hence a study of the in- 
variant geometrical properties of the black- hole merger, 
which would yield a highly-nontrivial perturbation of the 
'Kerr' background, may answer many open questions. 



II. MATHEMATICAL TECHNIQUES 

In the following sections we will use the convention 
that Latin indices range over the spatial coordinates (i.e. 
a = (1,2,3)) and Greek indices range over all four coor- 
dinates. 



A. Petrov type 

The Petrov classification of a generic spacetime is re- 
lated to the number of distinct principle null directions 
(PND) of the Weyl tensor. A generic spacetime will have 
four linearly independent null vectors fc^ (i.e. PNDs) at 
all points that satisfy 



k fc''fc[rC';^]iyp[o-^x] ~ ^' 



(2) 



Type I spacetimes have four distinct PNDs, Type II have 
three distinct PNDs (1 pair and two additional distinct 
PNDs), Type III have two distinct PNDs with one PND 
of multiplicity three. Type D spacetimes have two dis- 
tinct PNDs consisting of two pairs of PNDs of multiplic- 
ity two, type N spacetimes have a single PND of multi- 
plicity four, and Type O spacetimes have C^^pa = 0. 

If the tetrad is chosen such that is a PND, then the 
Weyl scalar -00 = C^,ypal^myiPmf vanishes, and simi- 
larly, if ijjQ = 0, then 1°^ is a PND. Hence the algebraic 
classification of the spacetime can be obtained by find- 
ing the number of distinct choices of for which -00 = 0. 
This amounts to finding the roots (and multiplicity of the 
roots) of the quartic equation (See Ref. [70], Eq. (9.5)) 



(3) 



where . . . , V'4 the Weyl scalars in an arbitrary 
tetrad, restricted only by the condition -04 7^ 0. This is 
equivalent to finding a tetrad rotation such that ■00 = Oi 
and if the root is repeated, then in this tetrad, V'l = 
(similarly if the multiplicity of the root is 3 or 4 then 
-02 — and ip'i — Q respectively). If, as in type D space- 
times, there are two pairs of repeated PNDs, then we 
can choose a tetrad where the only non-vanishing Weyl 
scalar is "02- It is important to note that the algebraic 
classification is done pointwise. A spacetime, as a whole, 



is of a particular type, if at every point the algebraic 
classification is of that type. 

In order to determine if the numerical spacetime is al- 
gebraically special (within the numerical errors) we fol- 
low J71,] a nd [70.1 , Ch. 4. We start by defining the scalar 
invariants |73 | 

/ = ^a/s^iC'"^^' and J = -ic„;3^5C75^.C''''"^. (4) 

where Cap-tS = j{Cap^s + ^^aPpvC'^^s) (i.e. 1/2 the 
conjugate of the self-dual part of the Weyl tensor Cap-ys)- 
If a spacetime has repeated principal null directions it 
is algebraically special. If this is the case, Eq. ([3]) has 
at least two repeated roots. In any case, Eq. ([3]) can be 
transformed into a depressed quartic (see Eq. ^ below) 
that, in turn, can be converted into a depressed nested 
cubic with roots y, which satisfy the condition 

y^-IV + 2J = 0. (5) 

Algebraic specialty then implies 

= 27 J^ (6) 

i.e. S* = 1 in Eq. [TJ For Types // and D the invariants 
/ and J are non-trivial, while for Types ///, iV, and O 
they vanish identically. 

For practical applications, it is convenient to write the 
invariants in terms of Weyl scalars in an arbitrary null 
tetrad 



/ = 37A2^ 

J = -i^l 



- 4?/>i'03 + 04-00, 

- V^o'04'02 + 2lpi1p3lp2 



Ml 



V'0'03- 



(7) 
(8) 



In order to completely determine the algebraic type we 
reduce Eq. ([3]), by changing to the variable a; = A'04 + 
"03 [ill, to the form 

x'^ + 6Lx^ + 4:Kx + N = 0, 



where 



K 
L 

N 



01 V'4 ~ 3V'4V'3V'2 
02-04 - ll'l, 



2^1 



3L^ 



04^Vo 



4V'4^V'lV'3 + 6V'402V'3^ 



(9) 

(10) 
(11) 

3^3^ (12) 



(note the typo in the definition of N in Refs. [70|, l7l|). 
For a type // spacetime, K and N — 9L^ ^ 0, while 
for type D and /// spacetimes, K = and TV — 9L^ = 
with iV 7^ 0. For a type N spacetime, K = and L — 
(hence = 0). 

Note that the above scalar objects are not invariant 
under arbitrary tetrad rotations (See Ref. (T^I, Chapter 
1, Eqs. (342), [note typo there], (346) and (347)). Tetrad 
rotations are classified as Type I, II, and HI, and have 
the form: 



(13) 
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+ on'^, 












-> An^", 






-> e'V^, 









which has the roots 



(14) 



(15) 



for Type I, II, and III, respectively, where a and b are 
complex scalars and A and are real scalars. Under 
these rotations the scalars L, K, and N transform as 



L - 


A^e 




K - 






N - 


A^e 





(16) 



for Type III rotations and 



L L, 
K K, 

N N, 



(17) 



for Type II rotations. Expressions for Type I rotations 
do not have these simple forms, but we verified that, if 
as in type D solutions, K = and N — 9L^ = in the 
original tetrad, then K ^ and N — 9L^ = in the 
new rotated tetrad (this is also obvious for type III and 
II transformations above). One the other hand L = is 
not preserved by type I rotations. 

Coming back to the roots xi,X2,X3,X4 of Eq. ([9]), we 
observe that, in numerically generated spacetimes, the 
roots never agree exactly, even if the metric is expected 
to be of a special algebraic type. Of course, the root dif- 
ferences in each pair should scale with resolution and 
asymptotically approach zero as h and t — s- oo 
(where h is the gridspacing). 

The roots of Eq. ^ can be obtained from the roots of 
Eq. ([H using the following algorithm [t^ 



D 
A 

yi 

2/2 
2/3 



J2- (7/3)3, 
A + B, 



B = {-J-Vd 



,1/3 



2^ 



B) 



1 



(A + B) 



^%A-B), 
^^{A-Bl 



(18) 



where the complex phases of A and B are chosen such 
that AB ^ 1/3. The roots of Eq. ^ are then obtained 
from the roots of the complete cubic equation for the 
variable z (where z — 2'04 y — AL) 

z^ + 12 Lz"^ + A{9L^ - N)z-IQK ^0, (19) 



zi = 2i>i yi - 4L, 
Z2 = 2V>4 y2 - 4L, 
Z3 = 2^-4 y3-4L. 



(20) 



Finally the roots of our original equation ([3]) can be writ- 
ten in the form [73>| 



Ai 

A2 

A3 
A4 



--03 + ^(-\/^+ \/z2 - V^) /V'4, 
-^3 + ^(-%Ar- V^ + Vi^) M, (21) 



where the signs of the ^/zi are chosen such that 
(\/^V^\/^) ~ ~4:K. We note that in a type D space- 
time Ai = A2 and A3 = A4. 



B. Vacuum 

The determination of the algebraic type of the matter 
fields can be done in an analogous way using the Ricci 
tensor, rather than the Weyl scalars. The analogue of 
the Petrov types are the Segre types and the equation 
to determine the multiplicities of the roots is ([TQi], Eq. 
(9.2)) 



1 



-her' 



ll7^ + l{ll-2h) = 0, (22) 



where 



and 



na qP C7 

Qa c/3 qS 
Dp a^j,. 



RaP — -gapR, 



(23) 
(24) 
(25) 



(26) 



is the trace free part of the Ricci tensor. 

This characterization of the matter fields does not com- 
pletely determine the algebraic properties, and other ad- 
ditional criteria have to be used. In our numerical simu- 
lations here, we are concerned with vacuum spacetimes. 
Numerical evolutions may introduce artificial (and un- 
physical) matter fields through violations of the Hamil- 
tonian and momentum constraints, and the natural way 
of monitoring the accuracy of the solution is to examine 
these constraints and confirm that the induced matter 
fields converge to zero. 



4 



C. Determination of the Kerr solution 

Once we determine that a solution is, for instance, 
Petrov type D and is a vacuum solution, we still do not 
uniquely single out the Kerr spacetime. One can go fur- 
ther and try to determine if the spacetime has the sym- 
metries of Kerr (the Kerr spacetime has two commuting 
spacelike and timelike Killing vectors [zBl)- However, one 
still needs to examine the asymptotic behavior of the so- 
lutions to determine that the spacetime does not have a 
NUT charge / or acceleration a. 

A general type D, vacuum Black hole solution can be 
described by the metric ([z3|, Eq- (17)), 



ds^ = i- 1-^ [di - (asin^ 6 + Al 
adt - (r'^ + (a + if'jdcf) 



sm 



P 

where 



2 o2 



and 



H. = 1 — a{l + acosO) r, 
= + {I + a cos 9)"^, 

P = sin^ 6(1- as cos 6 - a4 cos^ 6), 

Q = k — 2mr + er^ — 2anr^ — a^kr^ , 



as = 2aam~Aa^alk^ 
04 = —cP' 0? k 



(27) 



(28) 
(29) 
(30) 
(31) 



(32) 
(33) 



with e, 71 and fc as given a function of the more basic 
parameters to, Z, a, and a by 



+ Aalm - (a^ + 3P)a^fc, (34) 
■ - a[a^ -F)m+ {a^ - f )la^k, (35) 



kl 



-I' 



iaH^ fc = 1 + 2alm. 



(36) 



If the null tetrad is aligned with the principal null di- 
rections, i.e. 



)=(r^dr-da)+ VQdr 



(1 — apr) 
V2(r2+p2) 

(1 — apr) 
v/2(r2+p2) 

(1 — apr) 
" v/2(r2 +p2) 



then the only non-vanishing Wcyl scalar is 

3 



(37) 



^2^-{m + in) I L^PL \ . (38) 

r + ip 



where p = I + a cos 9. 

It is then natural to look at the asymptotic behavior 
of the spacetime to determine if there is a NUT charge Z, 
an acceleration a, or if the spacetime is plain Kerr. One 
can use the method of determining a quasi-Kinnersley 
frame jsoL Fts! ] to compute ^^2 and perform the above 
analysis. Alternatively, we can use the fact that, once 
we determined the spacetime is type D, we can choose a 
tetrad where all the Weyl scalars, but -02, vanish. Hence 
the invariants / and J must have the form 



(39) 



in this special class of tetrads. 

If the acceleration a 7^ then a series expansion of the 
invariant / gives 

/ = 3(TO+i/)^aV-— (m+i0^aV(«ap^ + l)+C ( 
r 

(40) 

Note that if the acceleration a = {) then n = I. An 
asymptotic expansion of the / invariant for the met- 
ric ([27l) then gives 

/ = -^(to -I- i/)^ - ^(m -I- il)'^{l + acosO) + O 

^41) 

and, by looking at the real and imaginary parts of the / 
invariant at large radii, we can determine the / parameter 
via 



2ml . 

—6a cos 



I ^ 0, 
I = 0. 



(42) 



We will use this method to determine the asymptotic be- 
havior of the final remnant of a black-hole-binary merger. 
Note that using / and J only requires smooth second 
derivatives of the metric, which has a distinct advan- 
tage over higher-derivative methods when dealing with 
numerically generated spacetimes. 



III. NUMERICAL TECHNIQUES 

To compute the numerical initial data, we use the 
puncture approach [7§| along with the TwoPuNC- 
TURES [131 code. In this approach the 3-metric on the 
initial slice has the form jab = {'4>bl + u)'^5ab, where 
ipBL is the Brill-Lindquist conformal factor, 6ab is the 
Euclidean metric, and u is (at least) on the punc- 
tures. The Brill-Lindquist conformal factor is given by 
il^BL = 1 + X]"=i "^^/(^l^— ^i|)j where n is the total num- 
ber of 'punctures', to^ is the mass parameter of puncture 
i {mf is not the horizon mass associated with puncture i), 
and fi is the coordinate location of puncture i. We evolve 
these black-hole-binary data-sets using the LazEv [81| 
implementation of the moving puncture approach 
In our version of the moving puncture approach we re- 
place the BSSN [H, m, [131 conformal exponent (j), which 
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has logarithmic singularities at the punctures, with the 
initially field x = exp(— 40). This new variable, along 
with the other BSSN variables, will remain finite pro- 
vided that one uses a suitable choice for the gauge. An 
alternative approach uses standard finite differencing of 
4> ■ Recently Marronetti et al. [85] proposed the use of 
W = y/x an evolution variable. For the runs presented 
here we use centered, eighth-order finite differencing in 
space jlll and a fourth-order Runge-Kutta time integra- 
tor (note that we do not upwind the advection terms). 

We use the Carpet [11] mesh refinement driver to pro- 
vide a 'moving boxes' style mesh refinement. In this ap- 
proach refined grids of fixed size are arranged about the 
coordinate centers of both holes. The Carpet code then 
moves these fine grids about the computational domain 
by following the trajectories of the two black holes. 

We obtain accurate, convergent waveforms and horizon 
parameters by evolving this system in conjunction with 
a modified l-|-log lapse and a modified Gamma-driver 
shift condition 0, l83|, and an initial lapse a{t — 0) = 
2/(1-1- The lapse and shift are evolved with 

{dt - P'di)a ^ -2aK, (43a) 
dtf}'' = B", (43b) 
dtB" = S/AdtT" -rjB". (43c) 

These gauge conditions require careful treatment of Xj 
the inverse of the three- metric conformal factor, near the 
puncture in order for the system to remain stable 
As shown in Ref. \8§\, this choice of gauge leads to a 
strongly hyperbolic evolution system provided that the 
shift does not become too large. In our tests, W showed 
better behavior at very early times {t < lOM) (i.e. did 
not require any special treatment near the punctures), 
but led to evolutions with larger truncation errors (im- 
portantly, larger orbital phase errors) when compared to 
X- 

We use AHFinderDirect [8^ to locate apparent 
horizons. We measure the magnitude of the horizon 
spin using the Isolated Horizon algorithm detailed in [6^ . 
This algorithm is based on finding an approximate rota- 
tional Killing vector (i.e. an approximate rotational sym- 
metry) on the horizon 93°. Given this approximate Killing 
vector if"", the spin magnitude is 

Sl<p]^^<f {^''R'Kab)d^V, (44) 
Stt J ah 

where Kab is the extrinsic curvature of the 3D-slice, (PV 
is the natural volume element intrinsic to the horizon, 
and i?° is the outward pointing unit vector normal to 
the horizon on the 3D-slice. We measure the direction of 
the spin by finding the coordinate line joining the poles 
of this Killing vector field using the technique introduced 
in 6]. Our algorithm for finding the poles of the Killing 
vector field has an accuracy of ~ 2° (see Q for details). 
Note that once we have the horizon spin, we can calculate 
the horizon mass via the Christodoulou formula (which 



is exact for a Kerr black-hole) 

= y'mf„ + 5V(4m2J, (45) 

where mi„ — yl/(167r) and A is the surface area of the 
horizon. 

We also use an alternative quasi-local measurement of 
the spin and linear momentum of the individual black 
holes in the binary that is based on the coordinate ro- 
tation and translation vectors [22] . In this approach the 
spin components of the horizon arc given by 

= '/'w^'^abdV, (46) 

J AH 

where (/ij^j = 6ij6rnkr™e'-^'^ , ^ ^ = 1, and r'" = x™ — 
is the coordinate displacement from the centroid of the 
hole, while the linear momentum is given by 

= ^ £^ (li]R''iKab - K^ab)d'V, (47) 

where ^j^j = 6}. 

A. Numerical Tetrad and Root Finder 

We calculate -00 ■ • • "04 using the tetrad 

r = + r>')/V2, (48) 
n^' = (t^-r^)/\/2, (49) 
= (61'' +i0'')/\/2, (50) 

where is the unit normal to the t = const slices and 
{r^, Qf^, are unit spacelike vectors (with time com- 

ponent equal to zero) constructed as follows [60]. We 
start with the unit vector 

= 0^, (51) 

where 0° = {—y,x,0} , v'^ = v"- / y/v°-v'^jab, and •jab is 
the spatial metric. We then find the unit vector in radial 
direction perpendicular to 0" 

r"=P, (52) 

where 

and r° = {x,y, z}. Finally, we obtain 

r = (54) 

where 

r = i^'^'ebcdct^'r^. (55) 

With this choice of tetrad ipQ ■ ■ ■ ip/^ are all non- vanishing 
for Kerr spacetimes when the specific spin a is non- 
vanishing. 
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TABLE I: Initial data parameters for the numerical evolution. 
The punctures have mass parameters mf, horizons masses 
mf, momenta ±p, spins Si, and the configuration has a total 
ADM mass A/adm = I.0000004M. 



m^/M 


0.37752 


mllM 


0.42452 


m? /M 


0.46298 


/M 


0.57872 


XI jM 


-0.75023 


X2/M 


0.58004 


yi/M 


1.11679 


j/a/M 


-0.89449 


zi/M 
Sf/M'^ 

Sf/M^ 


-0.16093 


Z2/M 

Sl/M^ 


0.20338 


-0.020765 


0.12106 


0.065806 


-0.05532 


0.054697 


0.16178 




-0.134735 


pK/M 


-0.21376 




-0.012323 







B. Initial Data 

To generate the initial data parameters, we used ran- 
dom values for the mass ratio and spins of the binary 
(the ranges for these parameters were chosen to make 
the evolution practical) . We then calculated approximate 
quasi-circular orbital parameters for a binary with these 
chosen parameters at an initial orbital separation of 50M 
and evolved using purely PN evolutions until the binary 
separation decreased to 2.3Af. The goal was to produce 
a binary that had no particular symmetries, so that we 
can draw general conclusions from the results, while also 
merging very quickly (within IbM of the start of the 
simulation), to reduce the computational expense. The 
initial binary configuration at r = 50M was chosen such 
that q = mi/m2 = 0.8, Si/ml = (-0.2,-0.14,0.32), 
and S2/ml = (-0.09, 0.48, 0.35)^his is the same basic 
configuration that we used in [90|- We summarize the 
initial data parameters in Table IH 

IV. RESULTS 

We ran the binary configuration using 9 levels of refine- 
ment with an outer grid of resolution h — 3.2M extend- 
ing to ±416M. The resolution on the finest grid was h = 
Af/80. We analyze the Weyl scalars in the region r < 5Af 
where we had a resolution oih < A//20. This calculation 
is non-trivial because the magnitudes of the Weyl scalars 
can be quite small (we need to analyze these scalars at 
very late times when the waveform amplitudes are quite 
small), requiring very-high overall simulation accuracy. 
We found that the isolated horizon formulae and the ra- 
diated energy and angular momentum both predict sim- 
ilar remnant masses and spins, with the isolated horizon 
formulae Eqs. dUl-dinD giving Mrem = 0.9859, S^om = 
{0.00160 ± 0.00005, 0.0407 ± 0.0004, 0.7173± 0.0001} and 
the radiation giving Mrcm = 0.9861 ± 0.0001, 5rcm = 
{0.00153 ± 0.00001, 0.04078 ± 0.00002, 0.7179 ± 0.0001}. 
A fit to the quasi-normal profile ~ exp(— af) sa\(Ljt) gives 
a = 0.07997 ± 0.0013 and lj = 0.5603 ± 0.0025, where 




100 200 300 



FIG. 1: The magnitude | (ys — 1/2) /yi \ versus time at the point 
X = 5M, 2/ = 0, 2; = 0. The spacetime is algebraically special 
if l(j/3 ~ y2)/yA ~ 0- Note the initial exponential decrease in 
the root difference. 



the values quoted are the average from fits to the real 
and imaginary parts of the (£ = 2,m = 2) compo- 
nent of ?/'4 extracted at r = lOOM over the domain 
(160M < t < 200M). The resulting values of M^cm 
and a/Mrem [iH are 0.9876 ± 0.0079 and 0.743 ± 0.013 
respectively. Note that the isolated horizon and radi- 
ated Energy/Momentum formulae predict that the final 
specific spin is a/M^em = 0.73931 ± 0.00016. This agree- 
ment is consistent with the final remnant being a Kerr 
hole (Note that this consistency is not a proof that the 
remnant is Kerr). 

If the spacetime is algebraically special, then the roots 
2/2 and 1/3 of Eq. ([5]) are equal. To measure how far 
the spacetime is from being algebraically special we plot 
the magnitude \(jj3 — y2)/yi\ (here yi provides a natu- 
ral normalization) and the invariant 5 — 1 at t he p oint 
{x = 5M,y = z = 0) (See Figs. [T] and [2]) [si, [salil. 

From the figures we can see that the deviation of the 
spacetime from being algebraically special decreases ex- 
ponentially (with an e-folding time of ~ 20Af for yi — 2/2 
and - lOAf for 5-1) with time until t - 150Af. The 
oscillation seen after this time may be due to refiections 
off of the refinement boundaries (this, in turn, provides 
a sensitive test to improve the numerical techniques). 

In Figs. [3][7] we show the unnormalized magnitudes of 
the root-pair differences |Ai — A2I and IA3 — A4I both as a 
function of i at a fixed (x, y, z) = (5, 0, 0) and along the x- 
axis at several times. Both pairs show a general decrease 
in the magnitudes of the differences with time, but with 
a pronounced oscillatory behavior. Note that |Ai — A2I 
separation is much smaller than the IA3 — A4I separation, 
indicating that the space-time first approaches Type II 
(and hence is algebraically special with — 1 ~ 0) be- 
fore settling to Type D. In Fig. [5] we plot the values of 
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100 200 300 

t/M 



FIG. 2: The magnitude 15* — 1| versus time at the point x = 
5M, 2/ = 0, 2 = 0. The spacetime is algebraically special if 
5=1. Note the initial exponential decrease in S — 1. 



the pairs (Ai,A2) and (A3,A4) on the complex plane at 
the point (5, 0, 0) for times t = 57, • • • , 166.25 in steps of 
0.59375. From the plots we can see how each of the two 
roots in the root pairs approach each other. In Fig. [5] we 
plot the magnitude of the root separations normalized 
by the difference between the average value of the roots 
in each pair (note that IA2 — A3I has an e-folding time 
of - 30M). It takes about 80M of evolution, or 65M 
post merger, until the larger normalized root separation 
falls below 1. Finally, in Fig. [6] we show the L2 norm 
of the root separations along the x and y axes restricted 
to 2M < \x\,\y\ < 5M and 2M < |y| < lOAf (the 
restriction to |y| > 2M is such that the black-hole 
interior is not included in the norm). The poorer con- 
vergence of the norm over the larger domain is due to 
numerical errors in the more coarsely resolved regions. 

In Fig. [5] we plot r\I\ versus r/M along the +?;-axis 
and along the line [x = Q^y = z). The leading-order 
term if a ^ and I — has an {apf" dependence, where 
p = I + acos9 (See Eq. ^7^). If ^ = then p = a cos 6*, 
and along the y-axis, p^ ^ 10^^ (the remnant spin is 
slightly misaligned with the z-axis), but along the line 
{x — 0,y — z), p^ ~ 0.028. From the data on the y-axis 
we can only conclude that al is very small. However, 
along the diagonal, p^ ~ .028 -I- 0.30Z, which provides 
evidence that both a and I are small. In Fig. [TO] we 
plot the function r3(/)/5R(/) versus M/x along the +x- 
axis for various times from t ~ 100 to i ~ 350M. It is 
clear from the plot that this function docs not tend to 00 
at larger r, which indicates that the NUT charge of the 
space time vanishes (i.e. given that we already found that 
a vanishes). Hence we can see good evidence that the 
spacetime is approaching Type D with zero NUT charge 
and zero acceleration, and hence is approaching a Kerr 
spacetime. 



le+00 




FIG. 3: The magnitude of the root-pair separation |Ai — A2 
versus time for the two roots close to A = at the point 
X = 5M, y = 0,z = 0. 




100 200 300 

t/M 



FIG. 4: The magnitude of the root-pair separation IA3 — A4 
versus time for the two roots furthest from A = at the point 
X = 5M, y = 0, z — 0. Note that there is no rapid decrease 
in the IA3 — A4I which indicates that the spacetime is not 
approaching type D as fast as it is approaching type II. 



We have confirmed that the constraints converge to 
zero for our code outside of the horizons. For this sim- 
ulation the constraint violations where of order 10~^ at 
the horizons, and dropped off steeply with radius. Con- 
vergence of the constraints is important to show that the 
spacetime remains a vacuum spacetime outside of the 
remnant horizons. 
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100 200 300 

t/M 



FIG. 5: The magnitude of the two root-pair separations 
normalized by the magnitude of the differences of the av- 
erage value of the roots in each pair |A(i_2) — -^{3,4)1, where 
A(i,2) = (Ai + A2)/2 and A(3,4) = (A3 + \4,)/2. 




100 200 300 400 

t/M 



FIG. 6: The L2 norm of the root separations versus time along 
the X and y axis for 2 < |a;|, \y\ < 5 and 2 < |a;|, \y\ < 10. The 
region containing the black hole itself was excluded from the 
norm. 



V. CONCLUSION 

We have provided a method to classify numerically 
generated spacetimes according to their algebraic prop- 
erties. This is based on the use of the coincidence of the 
principal null directions for algebraically special space- 
times. In particular, we focus on the final remnant of a 
generic-black-hole-binary merger, that, according to the 
'no hair' theorem, is expected to produce a Kerr black 




-6 -5 -4 -3 -2 -1 1 2 3 4 5 6 



x/M 

FIG. 7: The magnitude of the root-pair separation IA3 — A4[ 
along the a;-axis for several values of t. At first the root sepa- 
ration decreases significantly with time, but eventually stabi- 
lizes as numerical errors due to reflections off the refinement 
boundaries and other numerical sources of error begin to dom- 
inate. 

100 I ' , ' , ' 



25 



E 




u.u^ 

A * -0.02 -0.01 0.01 



A 

-100 I ' • ' • ' 

-100 100 200 

Re{X) 

FIG. 8: The locations on the complex plane of the roots 
Ai, ■ ■ ■ , A4 for f = 57, 57.59375, ■ ■ • , 166.25 at the point {x = 
5,y = 0, z — 0). The insets shows the last 107 points. Note 
that A4 has the largest scatter in time and that the separation 
of Ai and A2 is not distinguishable on the overall plot. Ini- 
tially, the points at different times are scattered, but converge 
to a fixed limit at late times. 

hole, and hence be of algebraic (Petrov) type D (i.e. that 
the four principal null directions agree in pairs) . We give 
a measure of the agreement by normalizing the numerical 
differences between two nearby roots of Eq. ^ with the 
average separation to the other root pair in the complex 
plane. 
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FIG. 9: r|7j as a function of rjM along the t/-axis and the 
diagonal line [x = = z). Note that the behavior indicates 
that r|/|— >Oasr— >oo, which indicates that the acceleration 
a vanishes. 
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1=178.125 
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CD 



0.5 



- y 



0.2 



0.4 0.6 

M/x 
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FIG. 10: The ratio r9(7) /»(/) as a function of M /x along the 
x-axAS. Note that the behavior indicates that 9(/)/5R(/) 
as r ^ oo (i.e. M/x — > 0), which indicates that the NUT 
charge vanishes. 



We have been able to verify this agreement to order 
10"'' and 10~^ for the two pairs respectively. We find 
that the agreement of the two roots in each pair improves 
with evolution time and only appears to be limited by un- 
physical boundary effects (from the refinement and outer 
boundaries). The late-time behavior of these two root 
pairs implies that the spacetime near the remnant first 
approaches an algebraically special type II (with one pair 
of roots and two distinct roots) and over longer timescales 
approaches type D. We also analyze the invariant asymp- 
totic behavior of the spacetime and do not find evidence 
for non-zero acceleration or NUT parameters. Thus, our 
simulations would suggest that the spacetime indeed ap- 
proaches Kerr, which incidentally, is also a strong test 
of the stability of the Kerr solution under large, generic 
perturbations within the timescales of the simulation. 

These results represent the first such tests for generic 
binary mergers using modest computational resources. 
This naturally suggests that further studies, perhaps also 
involving other numerical evolution methods, such as 
Pseudo-spectral [9l,|93| and multi-patch, multi-block [ol, 
l96j . be used to test the algebraic structure of the rem- 
nants of binary mergers. Finally, the algebraic structure 
of the remnants from the merger of more than two black 
holes (e.g. close-encounters [5^ [s^l of multiple black 
holes), while expected to have the same structure as the 
remnants of binaries, could conceivably have different al- 
gebraic structures. Thus it would be interesting to use 
these techniques to examine those remnants. 



Acknowledgments 



We thank S.Dain for insightful comments and M.Mars 
for bringing the C-metric to our attention. We grate- 
fully acknowledge NSF for financial support from grant 
PHY-0722315, PHY-0653303, PHY 0714388, and PHY 
0722703; and NASA for financial support from grants 
NASA 07-ATFP07-0158 and HST-AR-11763.01. Com- 
putational resources were provided by Lonestar cluster 
at TACC and by NewHorizons at RIT. 



[1] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005), gr- 
qc/0507014. 

[2] M. Campanelli, C. O. Lousto, P. Marronetti, and 
Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006), gr- 
qc/0511048. 

[3] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and 



J. van Meter, Phys. Rev. Lett. 96, 111102 (2006), gr- 
qc/0511103. 

[4] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 

Rev. D 74, 041501(R) (2006), gr-qc/0604012. 
[5] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 

Rev. D 74, 084023 (2006), astro-ph/0608275. 



10 



M. CampanoUi, C. O. Lousto, Y. Zlochowcr, B. Krish- 
nan, and D. Merritt, Phys. Rev. D75, 064030 (2007), [35 
gr-qc/0612076. 

M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 73, 061501(R) (2006). [36 
B. Brugmann et al., Phys. Rev. D77, 024027 (2008), gr- 
qc/0610128. [37 
J. G. Baker, J. Ccntrolla, D.-I. Choi, M. Koppitz, and 
J. van Meter, Phys. Rev. D 73, 104002 (2006), gr- [38 
qc/0602026. 

F. Pretorius, Class. Quant. Grav. 23, S529 (2006), gr- [39 
qc/0602115. 

F. Pretorius and D. Khurana, Class. Quant. Grav. 24, 

S83 (2007), gr-qc/0702084. [40 
J. G. Baker, J. R. van Meter, S. T. McWiUiarris, J. Cen- 
trella, and B. J. Kelly, Phys. Rev. Lett. 99, 181101 [41 
(2007), gr-qc/0612024. 

A. Buonanno, G. B. Cook, and F. Pretorius, Phys. Rev. [42 
D75, 124018 (2007), gr-qc/0610122. 

J. G. Baker et al., Phys. Rev. D75, 124024 (2007), gr- [43 

qc/0612117. 

M. A. Scheel et al., Phys. Rev. D74, 104006 (2006), gr- [44 
qc/0607056. 

J. G. Baker, M. Campanelli, F. Pretorius, and Y. Zlo- [45 
chower. Class. Quant. Grav. 24, S25 (2007), gr- 
qc/0701016. 

P. Marronetti et al.. Class. Quant. Grav. 24, S43 (2007), [46 
gr-qc/0701123. [47 
H. P. Pfciffcr et al., Class. Quant. Grav. 24, S59 (2007), [48^ 
gr-qc/0702106. 

U. Sperhake, V. Cardoso, F. Pretorius, E. Berti, and [49 
J. A. Gonzalez, Phys. Rev. Lett. 101, 161101 (2008), 
0806.1738. [50 
M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 
ritt, Astrophys. J. 659, L5 (2007), gr-qc/0701164. [51 
M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 
ritt, Phys. Rev. Lett. 98, 231102 (2007), gr-qc/0702133. [52 

B. Krishnan, C. O. Lousto, and Y. Zlochower, Phys. Rev. 
D76, 081501 (2007), 0707.0876. [53; 
S. Dain, C. O. Lousto, and Y. Zlochower, Phys. Rev. D 
78, 024039 (2008), 0803.0351. [54 
M. Campanelh, Class. Quant. Grav. 22, S387 (2005), 
astro-ph/0411744. [55 
F. Herrmann, D. Shoemaker, and P. Laguna, AIP Conf. [56 
873, 89 (2006), gr-qc/0601026. 

J. G. Baker et al., Astrophys. J. 653, L93 (2006), astro- [57 
ph/0603204. 

C. F. Sopuerta, N. Yunes, and P. Laguna, Phys. Rev. D [58 
74, 124010 (2006), astro-ph/0608600. [59 
J. A. Gonzalez, U. Sperhake, B. Brugmann, M. Hannam, 
and S. Husa, Phys. Rev. Lett. 98, 091101 (2007), gr- [60 

qc/0610154. 

C. F. Sopuerta, N. Yunes, and P. Laguna, Astrophys. J. [61 
656, L9 (2007), astro-ph/0611110. 

F. Herrmann, I. Hinder, D. Shoemaker, and P. Laguna, [62 
AIP Conf. Proc. 873, 89 (2006). 
F. Herrmann, I. Hinder, D. Shoemaker, and P. Laguna, [63 
Class. Quant. Grav. 24, S33 (2007). 
F. Herrmann, I. Hinder, D. Shoemaker, P. Laguna, [64 
and R. A. Matzner, Astrophys. J. 661, 430 (2007), gr- 
qc/0701143. [65 
M. Koppitz et al., Phys. Rev. Lett. 99, 041102 (2007), 
gr-qc/0701163. [66 

D. -I. Choi et al., Phys. Rev. D76, 104026 (2007), gr- 



qc/0702016. 

J. A. Gonzalez, M. D. Hannam, U. Sperhake, B. Brug- 
mann, and S. Husa, Phys. Rev. Lett. 98, 231101 (2007), 
gr-qc/0702052. 

J. G. Baker et al., Astrophys. J. 668, 1140 (2007), astro- 
ph/0702390. 

E. Berti et al., Phys. Rev. D76, 064034 (2007), gr- 
qc/0703053. 

W. Tichy and P. Marronetti, Phys. Rev. D76, 061502 

(2007) , gr-qc/0703075. 

F. Herrmann, I. Hinder, D. M. Shoemaker, P. Laguna, 
and R. A. Matzner, Phys. Rev. D76, 084032 (2007), 
0706.2541. 

B. Brugmann, J. A. Gonzalez, M. Hannam, S. Husa, and 
U. Sperhake, Phys. Rev. D77, 124047 (2008), 0707.0135. 
J. D. Schnittman et al., Phys. Rev. D77, 044031 (2008), 
0707.0301. 

K. Holley-Bockelmann, K. Gultekin, D. Shoemaker, and 
N. Yunes (2007), 0707.1334. 

D. PoUney et al., Phys. Rev. D76, 124002 (2007), 
0707.2559. 

I. H. Redmount and M. J. Rees, Comments on Astro- 
physics 14, 165 (1989). 

D. Merritt, M. Milosavljevic, M. Favata, S. A. Hughes, 
and D. E. Holz, Astrophys. J. 607, L9 (2004), astro- 
ph/0402057. 

A. Gualandris and D. Merritt (2007), 0708.0771. 

R. C. Kapoor, Pramana 7, 334 (1976). 

T. Bogdanovic, C. S. Reynolds, and M. C. Miller (2007), 

astro-ph/0703054. 

A. Loeb, Phys. Rev. Lett. 99, 041103 (2007), astro- 
ph/0703722. 

E. W. Bonning, G. A. Shields, and S. Salviander (2007), 
0705.4263. 

S. Komossa, H. Zhou, and H. Lu, Astrop. J. Letters 678, 
L81 (2008), 0804.4585. 

S. Komossa and D. Merritt, Astrophys. J. 683, L21 

(2008) , 0807.0223. 

G. A. Shields, E. W. Bonning, and S. Salviander (2008), 
0810.2563. 

L. RezzoUa et al., Astrophys. J. 674, L29 (2008), 
arXiv:0710.3345 [gr-qc]. 

U. Sperhake et al. (2007), arXiv:0710.3823 [gr-qc]. 

C. O. Lousto and Y. Zlochower, Phys. Rev. D77, 024034 

(2008), 0711.1165. 

M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 

Rev. D 77, 101501 (R) (2008), 0710.0879. 

P. O. Mazur (2000), hep-th/0101012. 

M. Campanelli, B. Kelly, and C. O. Lousto, Phys. Rev. 

D 73, 064005 (2006), gr-qc/0510122. 

,1. Baker, M. Campanelli, and C. O. Lousto, Phys. Rev. 

D 65, 044001 (2002), gr-qc/0104063. 

J. Baker and M. Campanelh, Phys. Rev. D 62, 127501 

(2000). 

O. Dreyer, B. Krishnan, D. Shoemaker, and E. Schnetter, 
Phys. Rev. D 67, 024018 (2003), gr-qc/0206008. 
M. Campanelli and C. O. Lousto, Phys. Rev. D 59, 
124022 (1999), gr-qc/9811019. 

C. O. Lousto and Y. Zlochower, Phys. Rev. D 76, 
041502(R) (2007), gr-qc/0703061. 

M. A. Scheel et al., Phys. Rev. D79, 024003 (2009), 
0810.1767. 

M. Mars, Class. Quant. Grav. 16, 2507 (1999), gr- 
qc/9904070. 



11 



[67] M. Mars, Class. Quant. Grav. 17, 3353 (2000), gr- 

qc/0004018. 

[68] B. F. Whiting, J. Math. Phys. 30, 1301 (1989). 

[69] G. Dotti, R. J. Gleiser, I. F. Ranea- Sandoval, and 

H. Vucetich (2008), 0805.4306. 
[70] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, 

and E. Herlt, Exact solutions to Einstein's field equations 

(Cambridge, UK: Univ. Pr., 2003), 2nd Edition. 701 P. 
[71] R. A. d'Inverno and R. A. Russel-Clark, J. Math. Phys. 

12, 1258 (1971). 
[72] J. Carminati and R. McLenaghan, J. Math. Phys. 32, 

3135 (1991). 

[73] L. Gunnarsen, H. Shinkai, and K. Maeda, Class. Quan- 
tum Grav. 12, 133 (1995), gr-qc/9406003. 

[74] S. Chandrasekhar, The Mathematical Theory of Black 
Holes (Oxford University Press, Oxford, England, 1983). 

[75] M. Abramowitz and I. Stegun, Handbook of Mathematical 
Functions (Dover, New York, 1972), 10th ed. 

[76] H. Stephani, Relativity : an introduction to special and 
general relativity (Cambridge Univ. Pr., Cambridge, UK, 
2004), 3rd ed. 

[77] J. B. Griffiths and J. Podolsky, Class. Quant. Grav. 22, 

3467 (2005), gr-qc/0507021. 
[78] C. Beetle, M. Bruni, L. M. Burko, and A. Nerozzi, Phys. 

Rev. D 72, 024013 (2005), gr-qc/0407012. 
[79] S. Brandt and B. Briigmann, Phys. Rev. Lett. 78, 3606 

(1997), gr-qc/9703066. 
[80] M. Ansorg, B. Briigmann, and W. Tichy, Phys. Rev. D 

70, 064011 (2004), gr-qc/0404056. 
[81] Y. Zlochower, J. G. Baker, M. Campanelli, and C. O. 

Lousto, Phys. Rev. D 72, 024021 (2005), gr-qc/0505055. 
[82] T. Nakamura, K. Oohara, and Y. Kojima, Prog. Theor. 



Phys. Suppl. 90. 1 (1987). 
[83] M. Shibata and T. Nakaumra, Phys. Rev. D 52, 5428 
(1995). 

[84] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 

024007 (1999), gr-qc/9810065. 
[85] P. Marronetti, W. Tichy, B. Brugmann, J. Gonzalez, and 

U. Sperhake, Phys. Rev. D77, 064010 (2008), 0709.2160. 
[86] E. Schnetter, S. H. Hawlcy, and 1. Hawke, Class. Quan- 
tum Grav. 21, 1465 (2004), gr-qc/0310042. 
[87] M. Alcubierre, B. Briigmann, P. Diener, M. Koppitz, 

D. Pollney, E. Seidel, and R. Takahashi, Phys. Rev. D 

67, 084023 (2003), gr-qc/0206072. 
[88] C. Gundlach and J. M. Martin-Garcia, Phys. Rev. D74, 

024016 (2006), gr-qc/0604035. 
[89] J. Thornburg, Class. Quantum Grav. 21, 743 (2004), gr- 

qc/0306056. 

[90] M. Campanelli, C. O. Lousto, H. Nakano, and Y. Zlo- 
chower (2008), 0808.0713. 

[91] F. Echevern'a, Phys. Rev. D 40, 3194 (1989). 

[92] J. Baker, M. Campanelli, C. O. Lousto, and R. Taka- 
hashi, Phys. Rev. D 65, 124012 (2002), astro- 
ph/0202469. 

[93] M. Boyle, L. Lindblom, H. Pfeiffer, M. Scheel, and L. E. 

Kidder, Phys. Rev. D75, 024006 (2007), gr-qc/0609047. 
[94] M. Tiglio, L. E. Kidder, and S. A. Teukolsky, Class. 

Quant. Grav. 25, 105022 (2008), 0712.2472. 
[95] E. Schnetter, P. Diener, E. N. Dorband, and M. Tiglio, 

Class. Quant. Grav. 23, S553 (2006), gr-qc/0602104. 
[96] B. Zink, E. Schnetter, and M. Tiglio, Phys. Rev. D77, 

103015 (2008), 0712.0353. 



